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ABSTRACT 

By regarding the vacuum as a perfect fluid with equation of state p = —p, de Sitter's cos- 

(N ! 

mological model is quantized. Our treatment differs from previous ones in that it endows the 

(N ■ 

vacuum with dynamical degress of freedom, following modern ideas that the cosmological term is 

00 : 

a manifestation of the vacuum energy. Instead of being postulated from the start, the cosmological 

O 

constant arises from the degrees of freedom of the vacuum regarded as a dynamical entity, and a 

C-t '. 

time variable can be naturally introduced. Taking the scale factor as the sole degree of freedom 
' of the gravitational field, stationary and wave-packet solutions to the Wheeler-DeWitt equation 
are found, whose properties are studied. It is found that states of the Universe with a definite 
value of the cosmological constant do not exist. For the wave packets investigated, quantum ef- 
fects are noticeable only for small values of the scale factor, a classical regime being attained at 
asymptotically large times. 

PACS numbers: 98.80.Hw , 04.60.Gw 
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1. INTRODUCTION 

Quantum cosmology is hopefully relevant to describe quantum gravitational effects in the very 
early Universe. In view of the nonexistence of a consistent quantum theory of gravity, minisuper- 
space quantization, which consists in "freezing out" all but a finite number of degrees of freedom 
of the gravitational field and its sources and quantizing the remaining ones, is expected to pro- 
vide general insights on what an acceptable quantum gravity should be like. This line of attack, 
initiated by DeWitt [1], has been extensively pursued to quantize model universes with different 
symmetries and varying matter content, and allows one to conceive theories of initial conditions 
for the wave function of the Universe [2]. Manifold schemes have been devised to quantize gravity 
coupled to matter in minisuperspace, the commonest of such quantization methods being those 
that rely on the Wheeler-DeWitt equation, advocate the quantization of only the conformal factor 
of the spacetime metric, or perform canonical quantization in the reduced phase space. 

In inflationary cosmology de Sitter's model plays a fundamental role, since it describes the 
phase of rapid expansion during which the vacuum energy dominates the energy density of the 
Universe, and gives rise to a term in the energy-momentum tensor that corresponds to a cosmologi 
cal constant. In modern cosmology the terms vacuum energy and cosmological constant are used 
almost synonymously [3]. It seems, therefore, of interest to study quantum aspects of de Sitter's 
cosmological model by treating the vacuum as a dynam ical entity. In such a treatment, the cos- 
mological constant should not be postulated from the start, but should emerge from the dynamical 
degrees of freedom of the vacuum. A possible way to achieve this is by regarding the vacuum as 
a perfect fluid with e quation of state p = —p. This approach appears to be fruitful, has several 
attractive features from the thermodynamic point of view, and leads to interesting consequences 
in inflationary cosmology [4]. The standard way of dealing with de Sitter's model in quantum 
cosmology [5] is highly questionable because it invol ves a system with a single degree of freedom 
and one constraint, so that, strictly speaking, the system has no degrees of freedom at all and 
is empty of physical content. The assignation of dynamical degrees of freedom to the vacuum 
circumvents this dif Acuity and renders our method distinctive in its ability to make room for the 
introduction of a time variable. 

Accordingly, we shall adopt Schutz's canonical formalism [6] which describes a relativistic fluid 
interacting with the gravitational field. This formalism is especially adequate for our purposes, 
inasmuch as it has the advantage of ascribing dynamical degrees of freedom to the fluid. As it will 
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be seen, Schutz's a ction principle is successful even in the case of the vacuum in the sense that 
the cosmological constant appears dynamically as a manifestation of the degrees of freedom of the 
fluid that acts as the vacuum. 

In the quantum realm the properties of de Sitter's model will be investigated on the basis of the 
associated Wheeler-DeWitt equation. Because the super-Hamiltonian constraint is linear in one of 
the momenta, the Wheeler-DeWitt equation can be reduced to a bona fide Schrodinger equation. 

This paper is organized as follows. In Section 2 a Hamiltonian treatment of de Sitter's model is 
developed on the basis of Schutz's canonical formalism, which is proved, in the case of the vacuum, 
to lead to the correct classical equations of motion. In S ection 3 the Wheeler-DeWitt equation is 
written down and is shown to take the form of a genuine Schrodinger equation for an appropriate 
form of the inner product. In order for the Hamiltonian operator to be self-adjoint its domain must 
be restricted to wave functions that obey certain boundary conditions. General sets of stationary 
solutions to the Wheeler-DeWitt equation obeying said boundary conditions are found. Then, in 
Section 4, normalized wave-packet solutions to the Wheeler-DeWitt equation ar e found, and their 
properties analyzed. Section 5 is dedicated to final comments. 

2. DYNAMICAL VACUUM IN DE SITTER'S COSMOLOGICAL MODEL 

The line element for a homogeneous and isotropic universe can be written in the Friedmann- 
Robertson- Walker form (we take c = 1) 



where <7y denotes the metric for a 3-space of constant curvature k = +1,0 or — 1, corresponding 
to spherical, flat or hyperbolic spacelike sections, respectively. 

The matter content will be taken to be a perfect fluid, and Schutz's canonical formulation of 
the dynamics of a relativistic fluid in interaction with the gravitational field will be employed [6] . 
The degrees of freedom ascribed to the fluid are five scalar potentials </?, a, (3, 6, S in terms of which 
the four-velocity of the fluid is written as 



ds 2 = g v \dx v dx 



A 



N{tfdt 2 + R(t) 2 a ij dx i dx j , 



(2.1) 



U v = - (w + aP,„ + 0S >U ) , 



(2.2) 
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where \x is the specific enthalpy. By means of the normalization condition 

g vX U v U x = -1 (2.3) 

one can express /i in terms of the velocity potentials. The action for the gravitational field plus 
perfect fluid is 

S=[ d 4 xy/=g~ + 2 [ d 3 xVhhijK ij + I d 4 xy/=g~p (2.4) 

JM JdM JM 

in units such that c = IQixG — 1. In the above equation p is the pressure of the fluid, is the 
scalar curvature derived from the spacetime metric g u \, is the 3-metric on the boundary dM of 
the 4-manifold M, and is the extrinsic curvature or second fundamental form of the boundary 
[7]. The surface term is necessary in the path-integral formulation of quantum gravity in order 
to rid the Einstein-Hilbert Lagrangian of second-order derivatives. Variations of the pressure are 
computed from the first law of thermodynamics. 

Compatibility with the homogeneous spacetime metric is guaranteed by taking all of the velocity 
potentials of the fluid as functions of t only. We shall take p — (7 — 1) p as equation of state for 
the fluid, where 7 is a constant and p is the fluid's energy density (we shall eventually put 7 = 0). 
In the geometry characterized by (2.1) the appropriate boundary condition for the action principle 
is to fix the initial and final hypersurfaces of constant time. The second fundamental form of the 
boundary becomes Kij = —hij/2N. As described in its full details in [8], after inserting the metric 
(2.1) into the action (2.4), using the equation of state, computing the canonical momenta and 
employing the constraint equations to eliminate the pair (0,pe), what remains is a reduced action 
in the Hamiltonian form 

S r = j dt{Rp R + <p Pip + S Ps - NU) (2.5) 

where an overall factor of the spatial integral of (deta) 1 ^ 2 has been discarded, since it has no effect 
on the equations of motion. The super-Hamiltonian 7i is given by 
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n = 




6kR + pZ iT 3(7_1) e s . 



(2.6) 



The lapse iV plays the role of a Lagrange multiplier, and upon its variation it is found that 
the super-Hamiltonian 7i vanishes. This is a constraint, revealing that the phase-space contains 
redundant canonical variables. 

For 7 = the super-Hamiltonian contains neither the fluid's degree of freedom ip nor its 
conjugate momentum p^, so that these canonical variables can be simply dropped. Equivalently, 
the correct classical equations of motion can be obtained without taking into account the degrees 
of freedom described by ip,a,f3 and 6, that is, they could have been disregarded from the start. 
It is a pleasant circumstance that only the physically meaningful entropy density S is relevant for 
7 = 0. The action reduces to 




(2.7) 



with 



n = 




- QkR + R 3 e s 



(2.8) 



This can be put in a more suggestive form by means of the canonical transformation 



T 




p T = 



(2.9) 



Then 




(2.10) 
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where 



n = 




-6kR + R 3 p T . 



(2.11) 



The extended phase-space is generated by (R,T,p R ,p T ). The variable T is such that the Poisson 
bracket 



so that T is a "global phase time" or, more precisely, since it does not involve the canonical 
momenta, a "global time" in accordance with the terminology introduced by Hajicek [9]. This is 
reassuring because the existence of a global time appears to be a necessary condition to prevent 
violations of unitarity in the quantum domain. 
The classical equations of motion are 



{T,H} H =o = R 3 >0 , 



(2.12) 



R = 



d(NH) 
dp R 



Np R 
12R 



(2.13a) 





T = 



d{NH) 
dp T 



= NR 3 , 



(2.13c) 



Pt = 



d{NH) 
dT 



(2.13d) 



supplemented by the super-Hamiltonian constraint 



7 



H = -- ^ - QkR + R 3 p T = . (2.14) 
24it 



In order to solve these equations in the case k — let us choose the gauge t = T, so that 
N = R~ 3 . Since px is constant, it follows from TC = that pr is proportional to R 2 . Insertion of 
this result in Eq.(2.13a) leads to 

R 2 R = constant =>• R(t) = (At) 1/3 , (2.15) 

where A is a positive constant and the origin of the time t has been conveniently chosen. The 
lapse function is, therefore, 

N(t) = iT 3 = (At)' 1 . (2.16) 



In terms of the cosmic time r defined by 



dr = N{t)dt = ^ =► ^-70 = ^ (2.17) 



one recovers the usual de Sitter solution 



R = R e Ar . (2.18) 



This concludes the verification that the action (2.10) leads to de Sitter's spacetime solution to 
Einstein's equations. Note that if the time variable is chosen as t = T the scale factor vanishes at 
t — 0, whereas in the cosmic-time gauge it vanishes only at r = — oo. It is also worth mentioning 
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that Eqs.(2.9) and (2.13d) show that the entropy density S remains constant, in agreement with 
the behavior of inflat ionary models during the de Sitter phase [3]. 

The general case of arbitrary k can be easily handled in the gauge N — 1 and leads to the 
expected solutions, but we shall refrain from considering it here. 

3. QUANTIZED MODEL: A WHEELER-DeWITT DESCRIPTION 

It will be convenient to introduce a new parametrization of the lapse function by writing it as 
NR. Then the action retains the form (2.10) but the super-Hamiltonian is now 

U = -^f - 6kR 2 + R 4 p T = . (3.1) 



The Wheeler-DeWitt quantization scheme consists in setting 



to form the operator H, and imposing the Wheeler-DeWitt equation 

H^ = (3.3) 

on the wave function of the universe In our present case this equation takes the form 

QkR ^ — iR = . (3.4) 



24 OR 2 dT 



Upon division by R A this equation takes the form of a Schrodinger equation 



1 d 2 ^ _ 6k 

l &f~UR^dR^'R^ ( ' 



with T playing the role of time. In order to be able to interpret T as a true time and (3.5) as a 
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genuine Schrodinger equation, the operator 



ft- 1 92 - 6 * (*R\ 

H ~2AR*dR? R 2 [6 - b) 



must be self-adjoint. The scale factor R is restricted to the domain R > 0, so that the minisuper- 
space quantization deals only with wave- functions defined on the half-line (0, oo). It is well-known 
that in such circumstances one has to impose boundary conditions on the allowed wave func- 
tions otherwise the relevant differential operators will not be self-adjoint. The need to impose 
boundary conditions to ensure self-adjointness has been long recognized by practitioners of the 
Arnowitt-Deser-Misner (ADM) reduced phase space formalism as applied to quantum cosmology 
[8,10-12], and very recently it has also been seen to have non-trivial cosmological implications in 
the Wheeler-DeWitt approach [13]. 

In the present case the operator H given by Eq.(3.6) with k = is self-adjoint in the inner 
product 

/•OO 

M>,<j>)= R 4 4>*(R)<p(R)dR (3.7) 
Jo 



if its domain is suitably specified. The operator H is symmetric if 

(i>, H<t>) = f rm^j^dR = jH ^j§^m)dR = m, 0) , (3.8) 

and, as in the case of d 2 /dR 2 on L 2 (0, oo), it is well known that the domain of self-adjointness of 
the Hamiltonian operator H comprises only those wave functions that obey 

^'(0) = o^(0) (3.9) 

with a G (—00, 00]. For the sake of simplicity, here we shall address ourselves in detail only to the 
cases a = oo and a — 0, that is, the boundary conditions we shall be mainly concerned with are 



*(0,T) = 



(3.10a) 
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or 



*'(0,T) = , 



(3.106) 



where the prime denotes partial derivative with respect to R. 



Let us look for stationary solutions to Eq.(3.4), that is, solutions of the form 



*(R,T) = e iET ij)(R) , 



(3.11) 



where E is a real parameter. Then the equation for ip{R) becomes 



1 d 2 j) 
24dR 2 ~ 



+ (ER 4 - 6kR 2 ) i) = . 



(3.12) 



The above equation coincides with the time-independent Wheeler-DeWitt equation written by 
other authors, occasionally with the help of somewhat obscure methods [14], with E playing here 
the role of the cosmological constant A. It should be emphasized that here this equation has 
been derived from a well-defined action principle and the cosmological constant has appeared 
dynamically from the vacuum degrees of freedom. 

In the de Sitter case (k = 0) it is easy to show from the above equation that the "cosmological 
constant" E is positive. Indeed, multiplying Eq.(3.12) by ip* and integrating over the half-line one 
finds 




(3.13) 



which, after an integration by parts followed by the use of (3.9), yields, for a > 0, 




(3.14) 



as we wished to prove. It should be clear from the above derivation that the general boundary 
condition (3.9) is not sufficient to allow us to reach the same conclusion. This special property of 
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conditions (3.9) with a > is not present in other minisuperspace models, and seems to confer 
this restricted set of boundary conditions a physically privileged status as compared to the general 
one with arbitrary a. 

The general solution to Eq.(3.12) with k = is [15] 



MR) = v r R\A J 1/e (PR 3 /3) + B J_ 1/6 (/?i? 3 /3) 



(3.15) 



where J v denotes a Bessel function of the first kind and order u, A and B are arbitrary constants, 
and 



(3 = V24E . 



(3.16) 



The usual interpretation of i? 4 |\l/| 2 as a probability density implies no correlation between R 
and T. The existence of such solutions to the Wheeler-DeWitt equation is perhaps not surprising 
since de Sitter's spacetime may be regarded as static [16] or self-similar [17]. 

It follows from the behavior of Bessel functions for small argument that in the case of boundary 
condition (3.10a) the solution is 

1>%\R) = VRJ 1/6 (PR 3 /3) , (3.17a) 

whereas in the case of boundary condition (3.10b) the solution is 

^ ] {R) = VRJ„ 1/6 ([3R 3 /3) . (3.176) 



From the asymptotic behavior of Bessel functions for small and large argument one easily checks 
that both solutions are square integrable, but their norm induced by the inner product (3.7) is 
infinite. Thus, states of the Universe with a definite value of the cosmological constant do not exist. 
Realizable states can only be constr ucted by superposition of solutions to the Wheeler-DeWitt 
equation with different values of the cosmological constant. 
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For any two states ip\ and ip2 belonging to the domain of the Hamiltonian operator, that is, 
obeying condition (3.9), one has 



Therefore, as in other minisuperspace models [12,18], here Vilenkin's wave function of the universe 
\l/ is ruled out because it is in conflict with the self-adjointness of the Hamiltonian operator. 
Indeed, Vilenkin's tunneling boundary condition [2,5] requires the wave function of the universe 
\1/ to consist only of outgoing modes at singular boundaries of superspace. In the present context 
this would amount mathematically to Ji2(0) > whenev er ipi = ip2 = which is impossible. 

4. EVOLUTION OF WAVE PACKETS 

The stationary solutions (3.17) have infinite norm and play here a role analogous to that of 
plane waves in the quantum mechanics of the free particle, that is, finite-norm solutions can be 
constructed by superposing them. The general solutions to the Whe eler-DeWitt equation (3.5) 
with k = are given by the continuous linear combinations 



where the superscript is used to distinguish the wave functions that obey the boundary condition 
(3.10a) from those that obey (3.10b). According to the Appendix, the probability distribution of 
values of the cosmological constant is given by 




(3.18) 




a = a,b 



(4.1) 



p^(E)= 1 -\c^(E)f , 



(4.2) 



assuming, of course, that ^^(R,T) is normalized in the inner product (3.7). 

We shall consider simple but illustrative examples of wave-packet solutions to the Wheeler- 
DeWitt equation obeying each of the boundary conditions (3.10). Introducing the parameter 
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we can write (4.1) as 



¥ a) (R, T) = VR r a W (A) e i3x2T/8 J„(\R 3 ) d\ 
Jo 



(4.4) 



where v — +1/6 or v — —1/6 according to whether a = a or b, and 



fl W(A) = ^cW(^) 



(4.5) 



The choice 



(4.6) 



with v = +1/6 for a = a and v = —1/6 for a = b, is particularly simple because it enables 
us to perform the integration in (4.4) and express the wave function of the Universe in terms of 
elementary functions [19]. In the case v — 1/6 we find 



*W(fl,D = [2(a - exp(- i ^-^ SIy ) , 



(4.7) 



whereas for z/ = —1/6 the result is 



r / 3iT, 
2(«-— ) 



-5/6 



exp (~4(^T)) 



(4.8) 



The expectation value of the scale factor is given by 
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/ °° R 4 ^*(R, T) RV(R, T)dR 
Jo 00 R 4 ty*(R,T)ty(R,T)dR 



(4.9) 



and its time dependence reflects the dynamical evolution of the quantized version of de Sitter's 
cosmological model. For the two types of boundary conditions of interest we find, respectively, 



(a) _r(|) r 64a 2 + 9TV/6 
\ n /T — 



32a 



(4.10) 



and 



r(|) 



64a 2 + 9T 2 



32a 



1/6 



(4.11) 



For sufficiently large values of the time T the expectation value (R)t grows at the same rate 
predicted by the classical solution (2.15), that is, the classical regime is attained for asymptotically 
large times. Quantum effe cts make themselves felt only for small enough T corresponding to small 
R, as expected. 

The dispersion of the wave packets defined by 



(ARf T 



(R 2 )t - (R) 2 T 



(4.12) 



is readily computed, with the results 



(AR) 



(a) 



/r(3/2) r(4/3) 



2\ 1/2 



\r(7/6) r(7/6) ; 



64a 2 + 9T 2 



32a 



1/6 



(4.13a) 



(AR)P = 



T(5/6)r(7/6) 

, r(5/6) 2 



1/2 



64a 2 + 9T 2 
32^ 



1/6 



(4.136) 
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The wave packets inevitably disperse as time passes, the minimum width being attained at T = 0. 
As in the case of the free particle, the more localized the initial state at T = the more rapidly 
the wave packet disperses. 

It is important to classify the nature of this model as concerns the presence or absence of 
singularities. For the states (4.7) and (4.8) the expectation value of R never vanishes, showing 
that these states are nonsingular. The issue of existence or nonexistence of singularities may be 
addressed from another point of view [20]. We can define the probability density 

P { °\R) = R A \^e\r)\ 2 , a = a,b , (4.14) 

for the stationary solutions (3.17a) and (3.17b). The behavior of the Bessel functions for small 
values of the argument makes it clear that p( a \R) — > as R — > 0, and thus the singularity 
is avoided within this mode 1 according to this criterion. Whatever the singularity criterion, de 
Sitter's quantum cosmological model is nonsingular just as its classical counterpart. 

5. CONCLUSION 

We have shown that taking the vacuum as a perfect fluid with equation of state p = —p a 
Hamiltonian description of de Sitter's cosmological model is possible, which makes subsequent 
quantization a straightforward process. This circumvents the problem of insufficient number of 
degrees of freedom that besets the usual Wheeler-DeWitt quantiz ation of de Sitter's model. The 
endowment of the vacuum with dynamical degress of freedom makes it possible the introduction 
of a time variable which, in turn, gives meaning to the dynamical evolution at the quantum level. 
The cosmological term is not postulated from the beginning, but arises as a manifestation of 
the vacuum degrees of freedom. In our approach states with a definite value of the cosmological 
constant are ruled out, and only those states are realizable that are finite- norm superpositions of 
solutions to the Wheeler-DeWitt equation with different values of the cosmological constant. 

With the scale factor as the sole degree of freedom of the gravitational field, stationary and 
simple wave-packet solutions to the Wheeler-DeWitt equation have been found. It turns out that, 
for the wave packets investigated, quantum effects are significant only for small values of the scale 



16 



factor, and the classical regime sets in at asymptotically large times. Just like the classical de 
Sitter model, its quantum counterpart is nonsingular. 
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APPENDIX: COSMOLOGICAL CONSTANT PROBABILITY DENSITY 

Let us start with a normalized state vector ty(R,t) for which 



roo 

m\ 2 = / R 4 \y(R,t)\' 

Jo 



dR 



roo roo roo 

/ dRR 4 / c{E)*e- lEt VRJ„(V24ER 3 /3)dE / c(E')e iE t VRJ u {V24E 1 R 3 /3)dE' . (Al) 
jo Jo Jo 



The change of variables 



E = ^\ 2 , E> = ^\' 2 , x = R 3 , g(X) = c(9A 2 /24) exp(-i9A 2 t/24) , (A2) 



leads to 



Q roo roo roo 

||^f = _ / dX\g(\y / dX'X'g(X') / xJ y {Xx)J v {X'x)dx . (A3) 
16 Jo jo jo 

With the help of Hankel's integral formula [21] 

roo roo 

f(x) = / J u (tx)tdt / f{X)J v (Xt)XdX , (A4) 
jo jo 

which is equivalent to the formal equation 
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/ xJJ\x)JJ\'x)dx = -5(X - A ) , (A5) 
Jo A 



one finds 

3 r 00 o 3 r°° 9A 2 , l0 1 f°° 



< roc -! roo MA f°° 

\m? = -l <^ W ? = -l d XX\c(-)? = -l \c(E)\ 2 dE , (A6) 



from which Eq.(4.2) follows. 
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